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We report a theoretical study of the thermopower of single-molecule junctions, with focus on 
phenyl-based molecular junctions. In contrast to prior studies, thermal fluctuations of the torsional 
angle between the phenyl rings and variations in the position of the molecular level alignment with 
respect to the electrode Fermi energy were taken into account. Full thermopower histograms were 
obtained, and their dependence on the magnitude of the fluctuations was studied. We found that 
at large molecular orbital variations the thermopower becomes strongly dependent on the torsion 
angle and can even change sign. This results in a marked effect of fluctuations on the thermopower 
distribution, yielding an average thermopower at high temperatures that differs (smaller or larger) 
from the fluctuation-free value, depending on the strength of fluctuations. We therefore conclude 
that fluctuations should be taken into account both when extracting single-molecule parameters, 
such as the molecular level-Fermi level offset, and in predictions of the thermopower of molecular 
junctions. 
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I. INTRODUCTION 

Current interest in the thermopower of single-molecule 
junctions is focused on two interrelated directions. The 
first is the expectation that molecular junctions may 
become the basis for high-efficiency thermoelectric de- 
vices, due to their low thermal conductivity and the 
large variety of possible junction compositions 1,2 . The 
second is the notion that the thermoelectric response 
carries inherent information regarding the mechanisms 
dominating the electronic transport and the electronic 
structure of the junction (primarily the positions of the 
molecule's highest occupied (HO) and lowest unoccu- 
pied (LU) molecular orbitals (MO) with respect to the 
electrodes' Fermi level) 3-5 . The relation between the 
two directions is simple - knowing the detailed elec- 
tronic properties of the junction will help design junc- 
tions with better thermoelectric performance. Against 
this background, the last few years have seen huge ad- 
vances in the measurement of thermopower in molecular 
junctions 6-12 , accompanied by a large number of theo- 
retical studies. 13-29 

The thermopower of a molecular junction (also known 
as the Seebeck coefficient) S is defined as the (linear) 
voltage response to an applied temperature difference, 
S = liniAT^o —AV/AT. Most of the theoretical papers 
mentioned above have used the non-equilibrium Green's 
function (NEGF) approach 30,31 in combination with den- 
sity functional theory (DFT) to determine S. Within this 
framework, the Kohn-Sham (KS) equations for the hy- 
brid metal-molecule-metal (finite) electrodes are solved 
self-consistently using various functionals 16 . The KS or- 
bitals are used to construct the Green's functions, which 
are, in turn, used to calculate the junction's transmission 
function, i.e., the probability of an electron arriving from 
the electrode with energy E to cross the junction. The 
self-energy that describes the electrodes is calculated ei- 



ther phenomenologically (using a simple level-broadening 
form 16 ) or by some self-consistent scheme 32 . Once the 
transmission function t(E) is calculated, it can be used 
to calculate the (temperature-dependent) thermopower 
S, as described in the following section. 

In recent years, the NEGF-DFT approach has been 
criticized, not only because of the approximation built 
into the use of KS orbitals to construct the Green's 
functions 31 or the use of different functionals 16 , but be- 
cause it cannot capture dynamical effects. Such effects 
may become very important for transport 33-35 , and their 
role in determining the thermopower is currently not 
known. In addition, it has recently been shown that be- 
cause of the inherent level broadening due to the self- 
energy, the NEGF approach "downplays" variations in 
the transmission function due to, e.g., the local den- 
sity of states at the molecule-electrode point of contact. 
As a result, fluctuations of the energy offset AE be- 
tween the HOMO and the electrodes' Fermi energy are 
overestimated 36 . 

The effect of fluctuations on the thermopower in molec- 
ular junctions has been a subject of recent theoretical 
studies 4,10,21,27,28,37-42 . However, neither the statistical 
nature (i.e. the full histogram of the thermopower) nor 
the effects of the variations in the molecular levels were 
discussed, even though all the experiments carry a signif- 
icant statistical signature. In fact, even in the much more 
thoroughly studied field of electron transport in molecu- 
lar junctions, the effect of various configurations on the 
statistical distribution were onlly recently discussed 43 . 
Addressing the effect of fluctuations theoretically is crit- 
ical for properly analyzing experimental results and re- 
lating them to the electronic properties of the molecular 
junctions. We thus conducted a theoretical study of the 
effect of fluctuations on the temperature-dependence of 
the thermopower in model single- molecule junctions com- 
posed of biphenyl (two rings) and triphenyl (three rings) 
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molecules attached to gold electrodes. 

Two sources of fluctuations were considered: (f ) fluc- 
tuations in the energy offset AE between the Fermi en- 
ergy of the electrodes and the HOMO level, and (2) ther- 
mal fluctuations of the torsion angle between the planes 
of the phenyl rings (the dihedral angle) 27 ' 44 . These 
two sources of fluctuations represent, respectively, fluc- 
tuations due to different reconstructions of the molec- 
ular junctions and thermal fluctuations within a given 
junction 8 ' 45 . Other sources of fluctuations, such as fluc- 
tuations in the molecule-electrode coupling may (and 
probably do) exist, and were discussed in past literature 2 . 
However, here we focus on the above sources of fluctua- 
tions, which have been experimentally proven to exist 45 
and to bare importance of the thermopower 8 . 

The main result reported here is that due to fluctu- 
ations, the average thermopower S(T) can deviate sub- 
stantially (higher or lower) from its bare, fluctuation free 
value. The full thermopower histogram is presented, and 
exhibits similarity to experimental results. The impor- 
tant conclusion that may be drawn from these findings 
is that in considerations of the thermopower of a molec- 
ular junction the full distribution of S(T) must be used 
to extract information on the molecular junction. 

II. MOLECULAR JUNCTIONS: MODEL AND 
CALCULATION 





eo t r 


S-TT1 


-4.0 -1.9 0.96 


NH2-TTI 


-4.30 -2.29 0.6 


CN-TT1 


-6.10 -2.0 0.14 



TABLE I: Parametrization of the TB model, following Burkle 



To address the above-described deficiencies, we ap- 
ply the 7r-orbital tight-binding (TB) model introduced in 
Ref. 50 and parametrized by Burkle et al. 27 . This model 
uses the transmission function from the full NEGF-DFT 
calculation as a fitting curve for a TB description of the 
biphenyl junction with different anchoring geometries. 
In the TB model, the molecule is described by a tt - 
orbital at every molecule location with an on-site energy 
eo and a hopping term t between orbitals belonging to 
the same phenyl ring. The inter-ring hopping is given by 
t' = t cos (p, where p is the torsion angle. In the wide- 
band approximation, the self-energy of the electrodes is 
characterized by a level broadening T/2. The TB model 
is schematically depicted in the top panel of Fig. 1, and 
the parameterization from Ref. 27 for different end- groups 
is shown in Table I. We focus here on the TT1 (top) 
binding geometry, as the differences between the differ- 
ent geometries are rather minor. 



In this work, we considered phenyl-based molecular 
junctions. Such junctions are prototypical molecular 
junctions whose conductance appears to be dominated 
by the torsion angle between the phenyl rings. 46-48 Note 
that the torsion angle can be adjusted by attaching to 
the rings different alkyl chains of different lengths 27 ' 47-49 . 
The equilibrium torsion angle po is bistable at po w 
30° and (po ~ 90°, with an energy barrier of about 
A « 0.3 eV between them (when the molecule is in 
the junction) 21 , 1 The thermopower of biphenyl molec- 
ular junctions was recently studied in detail with the 
NEGF-DFT method 21 ' 27 . Sergueev et al. 21 studied the 
effects of inelastic (electron-phonon) scattering on the 
thermopower for the two stable torsion angles. Although 
they did take into consideration thermal fluctuations of 
the torsion angle, they investigated only a single an- 
chor geometry between the molecule and the electrodes. 
Burkle et al. 27 considered different anchor geometries and 
different end-groups as well as different angles, but did 
not include the temperature dependence of the fluctua- 
tions in their study. Neither of these groups took into 
account the fluctuations of the HOMO level or the full 
thermopower distribution. 
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FIG. 1: Upper panel: Schematic representation of the molecu- 
lar junction and the tight-binding parametrization. (a) Trans- 
mission function t(E) (solid blue) and its logarithmic deriva- 
tive cr(E), which is proportional to the thermopower at low 
temperatures (dashed purple) as a function of energy E for a 
CN-TT1 configuration biphenyl junction. 



1 The energy barrier between the two stable states was calculated 
to be as be as small as 0.08 eV for free biphenyl molecules 44 , 
but we take here the higher limit as a " worse-case" . A smaller 
barrier would make the conclusions of this paper even stronger. 



Once the tight-binding Hamiltonian T~L{p) is 
parametrized (for a given torsion angle (p), we 
evaluate the retarded and advanced Green's 
functions G r,a (e, ip) = (e — H^irj) , from 
which the transmission function is calculated by 
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T(e,ip) = Tr[rG r (e)rG a (e)]. The thermopower can 
be calculated from r(e) by S{T) = -L x {T)/eTL Q {T), 
where L n {T) = J cLEt(E){E - [i) n [-df(E,T)/dE], 
e is the absolute value of the electron charge, 
f(E, T) = (exp[(E- n)/k B T] + ly 1 is the Fermi 
function, ks the Boltzmann constant, and /i is the 
chemical potential 2 ' 3 ' 31 . At low temperatures, the 
expression for S simplifies to S(E) = —q(T)a(Ep), 
where q(T) = n 2 k 2 B T/3e and a(E) = dlogr(E)/dE. 

In Fig. 1 we plot the transmission function r(e) (solid 
blue line) and the logarithmic derivative c(e) (dashed 
purple line) as a function of energy for the CN-TT1 con- 
figuration, which has the smallest level broadening. The 
Fermi energy for Au electrodes is taken at Ep ~ —5 
eV 27 (even though there is ambiguity of ±0.1 eV in 
this parameter 21 ' 29 ). The resonances at the HOMO and 
LUMO levels may be seen in the figure, with the reso- 
nance in a(E) close to them. 

In Fig. 2 (left panel) a 3D-map of the thermopower 
(at T = 40 K) is plotted as a function of energy eo, 
which determines the position of the HOMO level and 
the energy difference AE between the HOMO level and 
the Fermi energy (which is kept constant) and torsion 
angle tp (measured from the equilibrium value <po) . The 
thermopower is calculated using the full integral form and 
not with the low-temperature form. For e ~ —6.3 eV, 
we observe an increasingly strong dependence of S on the 
torsion angle. This strong dependence is evident from the 
upper panel of Fig. 2, showing a plot of S as a function 
of tp for two values of the energy offset, e = — 6.1eV 
(solid line, the value taken from Ref. 27 ) and eo = —6.5 
eV, for which the thermopower changes dramatically as 
a function of angle and even changes sign. The origin 
of this behavior can be understood from the observation 
that pushing eo to lower values is equivalent to raising the 
Fermi level from its bare value. From Fig. 1 we see that 
Ep = — 5 eV falls just below the resonance value at which 
the thermopower changes sign. Thus, pushing e down 
is equivalent to pushing the Fermi level closer to this 
resonance, where thermopower exhibits strong sensitivity 
to the torsion angle. 



III. FLUCTUATION-AVERAGED 
THERMOPOWER 

The strong sensitivity of S to the torsion angle begins 
to evolve within 0.2 eV of the bare value of eo- However, 
the value of eo (or, more physically, the energy offset AE 
between the HOMO with respect to the Fermi energy) 
is known to vary substantially between different recon- 
structions of the same junction 8 . Defining SE as the 
range of variations in AE (i.e., AE changes on a scale 
2SE between realizations) , we estimate SE to lie between 
a low of ~ 0.5 eV (extracted from transition voltage spec- 
troscopy measurement 45 ) and a high of SE ~ 1.4eV, as 
extracted from thermopower measurements 8 (although 
this may be an overestimate 36 ). Thus, taking into ac- 
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FIG. 2: Main panel: Thermopower S of the CN-TT1 con- 
figuration as a function of molecular orbital energy eo and 
biphenyl torsion angle tp. As eo drops below ~ —6.3 eV, 
a strong dependence of S on tp is observed. Upper panel: A 
comparison of S as a function of tp for two values of orbital en- 
ergy, eo = —6.1 and —6.5 eV, demonstrating the dependence 
of S on torsion angle. Note the possibility of sign change in 
S as a function of angle. 



count fluctuations of both eo and tp is essential in the 
calculation of the temperature dependence of the ther- 
mopower. 

To take both the variations in eo and the thermal fluc- 
tuations of tp into account, we calculate the average ther- 
mopower (S(T)) as a statistical average, taking eo from 
a uniform distribution U[e^ ± SE], where e^ are the 
bare values taken from Table I. It is reasonable to take 
a uniform distribution of eo, since the variations in e 
are due to junction reconstruction and are not thermal 
in origin. The phases are drawn from a thermal distri- 
bution T>(tp) oc exp ( — f~rji where £{tp) is taken phe- 

nomenologically to be a cosine function with two minima 
at tp — 30° and 90°, with an energy barrier AE = 0.3 
eV between them. Formally, the average thermopower is 
thus: 



(S(T)) 



$%-se de ° fo dL P s ( T > £ o, tp) exp (- 
2^/;d^exp(-g) 



k B T 



(1) 

In Fig. 3(a), we plot the average thermopower (S) as 
a function of the temperature, for the energy variation 
range SE = 0, 0.1, 0.6 eV (which are well within the 
experimentally relevant values 8,45 ), calculated for a CN- 
TT1 biphenyl molecular junction. The influence of the 
fluctuations can best be seen in Fig. 3(b), where the av- 
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FIG. 3: (a) Average thermopower (S) (averaged over tor- 
sion angles with a thermal distribution and molecular or- 
bitals with a uniform distribution) as a function of tem- 
perature T for different values of energy fluctuations range, 
SE = 0,0.1, 0.2, ...,0.6 eV. (b) (S) at T — 310 K as a func- 
tion of SE, demonstrating an inhomogeneous dependence, (c) 
Histograms of the thermopower S for SE = 0.1,0.3,0.6 eV. 
Note the double-peak structure (which was similarly found in 
experimental observations) and the negative values at large 
SE. 



erage thermopower (S) at T = 300 K is plotted as a 
function of the energy variation range SE. The average 
thermopower is not homogeneous in fluctuation strength, 
because as the fluctuations in eo become large, there are 
more frequent realizations for which eo is such that the 
Fermi energy is close to the transmission resonance, and 
thus the thermopower becomes strongly angle dependent 
(Fig. 2). At a certain fluctuation strength, the ther- 
mopower can even change sign with the torsion angle 
(for certain realizations), giving rise to negative ther- 
mopower values, which then tend to reduce the aver- 
age thermopower (hence the reduction at large SE and 
the non-monotonicity) . This phenomenon can be seen 
in Fig. 3(c), where the full thermopower S histogram is 
plotted for SE = 0.1, 0.3 and 0.6 eV. For SE = 0.6 eV, the 
histogram shows substantial weight on negative values of 
S. Note the histograms have a well-defined double-peak 
structure. Such structures, as well as sign-change of ther- 
mopower within the same junction, have been observed 
experimentally 6-8 . 

Out of all the configurations studied in Ref. 27 , in 
the CN-TT1 configuration of the biphenyl junction the 
molecule is most weakly coupled to the electrodes, as 
manifested by the smallest of all values of T (see Table I) . 
Thus, the CN-TT1 configuration has the sharpest trans- 
mission resonance and thermopower line-shape and, as a 
result, the strongest dependence of the thermopower on 
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FIG. 4: Averaged thermopower (S) as a function of orbital 
energy range SE at T = 310 K for the NH 2 -TT1 and S-TT1 
configurations. 



the torsion angle close to the resonance (Fig. 2). How- 
ever, the deviation of the average thermopower from the 
fluctuation-free value due to fluctuations in eo and (p is 
not limited to junctions exhibiting sharp resonance. In 
Fig. 4 we plot the average thermopower (S) at T = 310 K 
as a function of SE for the NH 2 -TT1 and S-TT1 config- 
uration biphenyl molecular junctions, with T = 0.6 and 
0.96 eV, respectively. Raising SE up to 1 eV (correspond- 
ing to a variation in the HOMO-fermi energy of 26 E = 2 
eV), we find for these junctions a strong dependence of 
(S) on 5E, including non-monotonic behavior and even 
a change of sign of (S) at large SE. 

Finally, we calculate the temperature dependence of 
a triphenyl molecule in the CN-TT1 configuration, de- 
picted in the upper panel of Fig. 5. The molecular junc- 
tion is now characterized by two torsion angles, ipi and 
tp2- To calculate the transmission and from it the ther- 
mopower, we use the same TB parametrization as that 
of the biphenyl junction, under the assumption that the 
hopping element, the orbital level and the coupling to 
the electrodes should not depend (or very weakly de- 
pend) on the number of phenyl rings. We allow the two 
torsion angles tpi and if2 to thermally fluctuate around 
their equilibrium values, and the average S(T, eo, <pi, ¥2) 
is averaged over the fluctuations of eo and the torsion 
angles. 

In Fig. 5, (S) is plotted as a function of temperature for 
an energy variation range of SE = 0, 0.2, 0.6 eV. Here 
we used the low-temperature form for the thermopower 
calculation. This is a good approximation, noting that 
the thermopower in Fig. 3 is almost linear in temper- 
ature, pointing that even room temperature is within 
the low-temperature regime (since the molecular energy 
scales are still much larger than the temperature). 

A monotonic reduction of (S) with SE is found, as 
seen in the inset, where we plot (S) as a function of SE 
at T = 300 K. Here we stress that for SE = 0.6 eV 
(a variation range of AE of 1.2 eV), the value of (S) 
is S = 17 /iV/K, which is ~ 25% of its value without 
fluctuations. If, for instance, we were to use this value of 
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FIG. 5: Average thermopower (S) as a function of tempera- 
ture T for different values of SE for the tri-phenyl molecular 
junction. Upper panel: Schematic representation of the tri- 
phenyl molecular junction with two torsion angles. Inset: Av- 
eraged thermopower (S) as a function of orbital energy range 
SE at T = 310 K. Note that for SE = 0.6 eV, the average 
thermopower is ~ 25% of its bare value, with no variations in 
e (i.e., SE = 0). 

S to calculate the position of the HOMO level without 
taking fluctuations into account, the resulting HOMO 
level would deviate by ~ 0.4 eV with respect to its real 
position. 

IV. SUMMARY 

In summary, we calculated the thermopower of phenyl- 
based molecular junctions. Using a DFT-based TB 
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